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ABSTRACT 

We present the results obtained by a detailed study of the extragalactic Z source, 
LMC X-2, using broad band Suzaku data and a large (~ 750 ksec) data set obtained 
with the proportional counter array (PCA) onboard RXTE. The PCA data allows 
for studying the complete spectral evolution along the horizontal, normal and flaring 
branches of the Z-track. Comparison with previous study show that the details of 
spectral evolution (like variation of Comptonizing electron temperature) , is similar to 
that of GX 17+2 but unlike that of Cyg X-2 and GX 349-f 2. This suggests that Z 
sources are heterogeneous group with perhaps LMC X-2 and GX 17+2 being member 
of a subclass. However non monotonic evolution of the Compton y-parameter seems to 
be generic to all sources. The broad band Suzaku data reveals that the additional soft 
component of the source modelled as a disk blackbody emission is strongly preferred 
over one where it is taken to be a blackbody spectrum. This component as well as 
the temperature of seed photons do not vary when source goes into a flaring mode 
and the entire variation can be ascribed to the Comptonizing cloud. The bolometric 
unabsorbed luminosity of the source is well constrained to be ~ 2.23 x 10"^^ ergs/sec 
which if the source is Eddington limited implies a neutron star mass of 1.6 M0. We 
discuss the implications of these results. 
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1 INTRODUCTION 

Low-mass X-ray binaries containing accreting weakly 
magnetised neutron stars exhibit luminosity and spec- 
tral variations on time scale of hours to days and are 
usually divided into two class es: Z and atoll sources 
(jHasinger and van der Klis|[l989h . This classification relies 
upon the path traced out by individual sources in the X-ray 
colour-colour diagram (CD) or hardness-intensity diagram 
(HID). CD is constructed by plotting the hard colours (ra- 
tio of count rates in two largest energy bands) against soft 
colours (ratio of count rates in two lowest energy bands). 
HID is constructed by plotting hard colour against inten- 
sity. Z-sources trace out an approximate 'Z' shaped path 
in the CD and HID with the three branches named as the 
horizontal branch (HB), the normal branch (NB) and the 
flaring branch (FB). Spectral analysis of UV data which 
is modelled as reproduced X-ray emission from an outer 
disk, indicates that the bolometric luminosity and hence the 
mass accretion ra te increases along the Z-track from HB t o 
FB through NB (jHasinger et al.lll990l : IVrtilek et al.lll99d) . 
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However, more direct estimation of the bolometric luminos- 
ity has proved to be difficult since that requires both broad 
band data and a reliable distance estimate. Nevertheless, 
there are indications that Z-sources have luminosities com- 
parable to the Eddington limit. There are nine Z-sources 
known, including the recent discovery of a transient galac- 



tic X-ray binary dHoman et al. 


20071). Two of them LMC 


X-2 JSmale and Kuulkerd 200C : 


Smale. Homan & Kuulkerd 


l2003f) and RX J0042.6+4115 (1 


3arnard et al.ll2003ll are ex- 



tragalactic ones. 



The radiative processes that produce the X-ray spec- 
tra of neutron star LMXBs are not well constrained. While, 
there have been attempts to have detailed theoretical model 
based on the accretion flow (e.g. [Psaltis Lamb & Milled 
Il995h . in general, two different phenomenological approaches 
have been adopted to model their spectra. In both these 
approaches, the spectrum is modelled as sum of two main 
components, one arising from an accretion disk and the other 
from the boundary layer between the disk and the neutron 
star surface. In the first approach (sometimes referred to 
as the "Western" approach), the accretion disk is consid- 
ered to be a standard cold disk emitting a soft component, 
while the boundary layer radiates the harder Comptonized 
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Table 1. RXTE observations of LMC X-2 



OBSID 




Start Date 


End Date 


Duration 










(ksec) 


50041-01-01- 




2001 Fob 10 


2001 Feb 14 


145 


60027-01-01- 




2001 Aug 30 


2001 Sep 3 


144 


60027-01-02- 




2001 Dec 13 


2001 Dec 16 


122 


60027-01-03- 




2002 Feb 01 


2002 Feb 06 


148 


70017-01-**- 




2002 Aug 22 


2002 Oct 30 


185 



component (iMitsuda et"apll984 Il989l : iDi Salvo et "al~ll2002l : 
lAgrawal and Sreekumai] 20031 ). In the second approach 
(called the "Eastern" approach), the boundary layer emits 
a cool black body emission while the h arder Comptonized 
component ar ises from an hot inner disk l|White et al. |[l98^ : 
iDi Salvo et al. 2000 . 2001i l . There is also evidence for a vari- 
able har d power-law component in the spectra of five Z- 
sources (lAsai et al.lll994l : lOi Salvo et"ani200d . l200ll . I2OO2I : 
iD'Amico et al. I 200ll ). The two approaches are difficult to 
distinguish primarily due to the absence of good quality low 
energy (0.1-1.0 keV) data and the uncertainty in the absorp- 
tion column density. 

To understand the phenomenology of these sources, it 
is important to study their complete spectral evolution in 
order to confirm any generic behaviour and to see if there 
are sub-classes of sources that have similar spectral evolu- 
tion. There is already some evidence that Z-sources may 
no t be an ho mogeneous group. T his was first indicated 
by I KuulkerTfc van dor Klis (1 99^ : 1 Kuulkers et all (119971 ) 
who classified them into 'Sco' and 'Cyg' like sources. De- 
tailed study of spectral evol ution along the Z-track has bee n 
carried o ut for GX 349-1- 2 ( Agrawal and Sreekumaij 20031), 
Cyg X-2 (|Di Salvo et al . "2002"), GX 17+2 ('Pi Salvo et al ' 



200c ) and GX 340-1-0 ([Church, Halai fc Bl aucinska-Church 
20061 ). In general, there is considerable variation of the spec- 
tral parameters as a source moves along the Z-track and the 
variations are not always similar in all sources, indicating 
that Z sources are not an heterogeneous group. For exam- 
ple, as Cyg X-2 moves from the HB to NB, the optical depth 
of the Comptonized component decreases while the elec- 
tron temperature increases dramatically from 3 to 8 keV. 
A more modest increase in electron temperature from 2.7 to 
3.2 is seen in GX 3494-2. On the other hand, for GX 17-1-2 
the temperature decreases along the same track. Along the 
fiaring branch, detailed spectral e volution study has only 
been undertaken for GX 349-1-2 (|Agrawal and Sreekumaj 
|2003| ) where the temperature increases along the track. 
While these results indicate that the physical parameters 
and conditions for these sources are different, what seems 
to be generic is that the effect of Comptonization decreases 
along the normal branch while it increases along the flaring 
branch. This is revealed by the variation of the Comptoniza- 
tion parameter y = A{kT,,/me(?)T^, where kTe is electron 
temperature, me is rest mass of electron, c is speed of light, 
and r is optical depth of corona. The Comptonization pa- 
rameter y decreases along the normal branch and increases 
along the flaring one, a behaviour clearly seen in GX 349-1-2. 

The temporal properties of Z-sources are also known 
to correlate with the position of the source in the Z-track. 
Three different kinds of quasi-periodic oscillations (QPOs) 
have been seen in these sources. QPO with frequency 15-100 



Hz, which appear in the horizontal branch (but sometimes 
seen up to the upper part of the normal branch), are called 
horizontal branch oscillations (HBO). The frequency of an 
HBO increases from top left to bottom right of the hor i- 
zontal branch (jHoman et al. II2002I : IWiinands et al. Ill997l '). 
QPO with frequency 5-15 Hz are observed in the normal 
and the flari ng branch, and are c alled normal-flaring branch 
oscillations ([Homan et al. II2OO2I ). These low frequency os- 
cillations have been seen in all the galactic Z-sources ex- 
cept G X 349-1-2 where instead broad pea k ed noise is ob- 
serve d (jAgrawal and Bhattacharvval l2003l : I O'Neill et al. I 
|2002| ). The third type of QPO are kHz ones (300-1200 
Hz) which are observed in all the branches. The frequency 
of kHz QPO increa ses along the Z-track from HB to FB 
(| van der Khsll2000D 

Early observation of LMC X-2 show ed that its lumi- 
nosit y varies from 0.6 to 3x 10^® erg s"^ (iMarket fc Clark I 
Il975h . LMC X-2 was identified as a Z-source when it traced 
out a complete Z-track during 2001-2(3 02 RXTE observa- 
tions (|Smale. Homan fc Kuulkerd 120031). Its optical coun- 
terpart is a faint blue star ( Pakull I 19781 ) typical of other 
LMXB. Simultaneous optical and X-ray observations of this 
source have revealed that X-ray and opt ical emission are 
corre lated with a time delay of ~ 20 s (|McGowan et al.~l 
|2003| ). A brief (~ 10 ksec) XMM-Newton observation of this 
source can be described by a disk blackbody, with inner disk 
temperature Ti„ ~ 0.5 keV and a blackbody componen t 
with temperature T^b ~ 1.5 keV (|Lavagetto et al. I I2OO8I ). 
However, the source is bright enough to cause pile up in 
XMM-Newton observations, making such spectral analysis 
difficult. No QPO has been detected in LMC X-2, perhaps 
due to low statistics data. 

LMC X-2 provides an unique opportunity to under- 
stand the spectral evolution of Z sources and the radia- 
tive processes that occur in such sources. The distance 
(^ 50 ± 2 kpc) to the sourc e (|Freedman et al. I I2OOII : 
iMcChntock and RemiUardI |2003| ) is much less uncertain 
than those of galactic sources enabling a more accurate esti- 
mation of its intrinsic luminosity. The Galactic column den- 
sity along the line of the source is 6.3 x 10^° cm~^. This 
provides a firm lower limit to the absorption which can pro- 
vide constrains on spectral models. The source has been ex- 
tensively observed by the Rossi X-ray Timing Experiment 
(RXTE) for a good time duration of ~ 750 ksec. This al- 
lows the construction of elaborate CCD and HID diagrams 
and a detailed study of its spectral evolution. In addition, 
a ~ 50 ksec observation by Suzaku reveals the broad band 
spectrum (0.3-30 keV) of this source which allows models to 
be differentiated and the bolometric luminosity to be esti- 
mated. 

In this paper, we report the results of spectral analysis 
of RXTE in §2 and Suzaku data in §3. In §4, we summarise 
the main results and discuss the implications. 



2 RXTE OBSERVATIONS 

2.1 Colour-Colour and Hardness Intensity 
Diagrams 

LMC X-2 was observed by RXTE during 2001-2002 for 
a total good time of 750 ksec. The Observation log is 
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Soft Color 

Figure 1. X-ray colour-colour diagram obtained using the RXTE-PCA data for five observations a) 2001 Feb 10 b) 2001 Aug 30 c) 2001 
Dec 13 d) 2002 Feb 01 e) 2002 Aug 22. Soft colour is ratio of count rates in the 4.5-6.5 and 2.5-4.5 keV energy bands and hard colour 
is ratio of count rates in the 9.8-18.5 and 6.5-9.8 keV bands. The colour values have been calculated using 512 seconds averages. The 
boxes mark the regions for which spectral analysis was undertaken. The representative error-bar for colour values is given in upper right 
corner of each figures 



shown in Table [T] We have analysed the data collected 
by Proportional-Counter-Array (PCA) which operates in 
2-60 keV energy band and consists of five proportional 
counter units with a total effective area of ~ 6500 cm'^ 
l|jahoda et al. Ill996l '). The analysis was done using PCA 
standard-2 mode data. Since only two of the five propor- 
tional counter units (PCUO and PCU2) were reliably on 
all through the observation, only data from these PCUs 
were selected for analysis. We define soft colour as ratio 
of count rates in the 4.5-6.5 keV and 2.5-4.5 keV energy 
bands and hard colour as that in the 9.8-18.5 keV and 



6.5-9.8 keV bands. The CD, generated using 512 sec av- 
erage counts, for different observations are plotted in Fig- 
ure [T] In Figure [5J a), the combined CD for all the obser- 
vations is plotted. The resu l t is si milar to that obtained by 
'Sma le. Homan fc Kuulkerj (120031 ) . except that the present 
analysis is for a larger data set. We define overall colour 
as ratio of count rates in the 6.5-18.5 keV and 2.5-6.5 keV 
bands. The Hardness Intensity Diagram (HID) was con- 
structed by plotting the overall colour against the 2.5-18.5 
keV counts/sec/2PCUs and is plotted in Figure [2jb). Since 
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Figure 2. a) A single, combined colour-colour diagram for the five different observations (see Table [TJ. The energy bands and time bin 
size used to construct the combined colour-colour diagram is same as that used in Fig. [T] b) Combined hardness-intensity diagram for 
five different observations. The overall colour is the ratio of count rates in the 6.5-18.5 and 2.5-6.5 keV energy ranges. The intensity is 
for 2.5-18.5 keV and is in units of counts/sec/2PCUs. Each points corresponds to 512 sec averages. The representative error bar for each 
data points is given in the upper right corner of each figure. 



there was no gain changes after 2000 there was no need to 
normalise the intensities and colours. 

As seen in Fig[2ja), the horizontal branch (HB) for this 
source is neither horizontal nor distinct. The classification 
requires study of other spectral relationships like HID which 
is shown in Fig[2l;b). The HID of LMC X-2 is quite simi- 
lar to that s een for another well studied Z-source GX 17-f 2 
(see Fig.l of lHoman et al. Il2002l l. In the HID of GX 174-2, 
the slanted branch (almost vertical) at the top left is called 
the horizontal branch and we do a similar identification for 
LMC X-2. We divide the HB in the HID, into two region 
HBl and HB2 and mark the corresponding location in the 
CD. The normal and flaring branch are distinct in the CD. 
Thus, we are able to identify all three branches for LMC 
X-2. For spectral studies, the Z-curve has been divided into 
eight regions. Two for the horizontal branch (HBl and HB2), 
two for the normal branch (NBl and NB2), a normal/flaring 
branch vortex point (N/FBV) and three regions for the flar- 
ing branch (FBI, FB2 and FB3). All eight regions have been 
marked out in Fig[2ja). 



2.2 Spectral Evolution 

We extracted spectra corresponding for the eight regions in 
the CD to carry out a a detailed spectral evolution study. 
We added 1 % systematic error to all the spectra to take in 
to account the uncertainty in the PCA response matrix. 

A single c omponent Compt onization model (compTT 
in XSPEC; see iTitarchuk |[l993 l with absorption provided 
adequate fit to all the PCA spectra in the energy range 3-20 
keV. Since the data is for energies > 3 keV, the absorbing 
column density is not well constrained. Hence we adopt a 
constant value of N_r/ 9 x 10^° based on the Suzaku analysis 
reported in section 3. An iron emission line centred at ~ 
6.4 is also required to improve the fit. In the top panel of 
Fig[3]we show the observed count rate spectrum for HBl, 
along with the best fit model. The residuals (in units of a), 
with and without inclusion of systematic error are shown in 
the middle and bottom panels respectively. The PCA data 



suggests that the soft component, originating either from 
the disk or from the boundary layer, is absent above 3 keV. 

The spectral parameters of compTT component evolve 
significantly as the source moves along the Z-track. The best 
fit parameters at different parts of Z-track have been listed 
in Table [2] and evolution along the track has been plotted 
in Fig. |4] The electron temperature decreases as the source 
evolves from HB to lower part of FB . This variation is differ- 
ent f rom that seen in GX 349-1-2 (jAgrawal and Sreekumaij 
l2003h and Cyg X-2 (|Di Salvo et al. I I2OO2I ) where the op- 
posite behaviour was found. On the other hand it is simi- 
lar to that observed in GX 17-1-2 (|Di Salvo et al. II2OO0I '). As 
LMC X-2 moves along the flaring branch the temperature in- 
creases which again is in contrast to the behaviour observed 
for GX 349+2 where the temperature decreases. However 
what seems to generic is that effect of Comptonization de- 
creases along the normal branch while it increases along the 
flaring branch. This is revealed by variation in Comptoniza- 
tion parameter y=4(kTe/meC'^)T'^ which decreases along the 
normal branch and increases along the flaring one (Fig. Q a 
behaviour also seen in GX 349-1-2. In the flaring branch as y 
increases the spectrum is that of saturated Comptonization 
with peak ~ 3kTe ~ 6 keV. 



3 SUZAKU OBSERVATION 

3.1 Colour-Colour Diagram and Light curves 

Suzaku observed LMC X-2 for a good time duration 
of 119 ksec on 2006 A pril 24 and 25. The observa- 
tory (Mitsuda et al.l |2007| ) consists of two different types 
of instrument , the X-ray Imaging Spectrometer (XIS; 
[Kovama et~aLl l2007h and Hard X-ray Detector (HXD; 
T akahashi et al.ll2007l ). 

There are four XIS units (named XISQ, XISl, XIS2, 
XIS3), each of them is a 1024 x 1024 pixel CCD and covers 
an energy range of 0.3-12 keV. HXD consists of two types of 
detectors, the Si diode detector (PIN) and the Gadolinium 
silicate crystal (GSO). The PIN diode covers an energy range 
of 10-70 keV while the GSO scintillators are sensitive in 40- 
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Table 2. The best fit parameters obtained by fitting the RXTE-PCA data with XSPEC model compTT. The electron temperature kTe 
and seed photon temperature kTw is measured in keV. r is optical depth of central corona and y is Compton y-parameter. The division 
N/FBV means normal-fiaring branch vertex 
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Figure 3. The observed PCA spectrum of LMC X-2 in the HB 
together with the best fit diskbb+compTT model are shown in 
panel 1. Residuals in unit of sigma without and with inclusion of 
systematic errors are shown in panel 2 and panel 3 respectively. 



Figure 4. The variation in spectral parameters obtained by fit- 
ting the PCA data with Comptonization model. The main vari- 
ation is that of the temperature of Comptonization cloud which 
decreases from the upper part of the HB to the lower part of the 
NB. The Compton y-parameter also show a similar variation 



600 keV band. The XIS data were collected in normal mode 
with 1/4 window option. The minimum bin time available 
in this mode is 2 seconds. The cleaned XIS event files were 
used to obtain the light-curves and spectra. The XIS image 
is shown in Fig. [5] The source light-curves were extracted 
using 4' radius centred at source position. 

To compare with RXTE observations, we study the 
colour and intensity properties of the counts at energies 
> 3 keV. Fig. [6] shows the colour-colour diagram for the 
Suzaku data using three energy bands defined as 3-4.5 keV, 
4.5-6.5 keV and 6.5-10 keV. The data points corresponds 
to 512 sec averages. Superimposed on the figure is the simu- 
lated colour-colour diagram, obtained using the best fit PCA 
spectral parameters and XISO response matrix. In particu- 
lar, eight fake X-ray spectra were created using the best fit 
RXTE model parameters (Table ^ using the Suzaku XIS 
response matrix. The count rates in three energy bands (3- 
4.5, 4.5-6.5, 6.5-10 keV) were calculated for these eight sim- 
ulated spectra and the simulated colour-colour diagram was 
computed. The simulated CD is represented by a solid line 
connecting the eight points. While it is difficult to differenti- 
ate the different branches, the Suzaku data points are consis- 




Figure 5. The XIS image of LMC X-2 
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Figure 6. The colour-colour diagram constructed using XISO 
data. The hard colour is ratio of count rates in 6.5-10 keV and 
4.5-6.5 keV energy bands, while the soft colour is ratio of count 
rates in the energy bands 4.5-6.5 keV and 3.0-4.5 keV. The solid 
line connects the eight simulated colours on the parameters cor- 
responding to RXTE observations (Table [21l. 



tent with simulated ones. Note that statistical error for these 
data points are significantly larger than for RXTE data be- 
cause the 3-10 keV count rates for Suzaku (~ 6 counts/sec) is 
smaller than that for RXTE, - 60-150 counts/sec/2PCUs. 
Comparison with hardness intensity plot generated using 
RXTE data is more difficult because apart from the statis- 
tical uncertainty of Suzaku points, the absolute calibration 
between the PCA and Suzaku may differ. 

The XISO fight-curve created using 3-10 keV energy 
band is shown in the top panel of Fig. [T] There are clear 
flaring like activity where the rate increases by ~ 30 %. 
Moreover the spectra hardens during the flares as revealed 
in the bottom panel of Fig. [7] where overall colour is plot- 
ted as function of time. The overall colour is ratio of count 
rates in 5-10 and 3-5 keV energy bands. Thus, we divided 
the data in two states, a flaring state where the count rate is 
> 7 counts/sec and a non flaring one corresponding to rates 
< 7 counts/s. The low statistics of the data does not allow 
for more reliable divisions. We identify the flaring state with 
the flaring branch and the low count rate state as the normal 
branch. 



3.2 Spectral Evolution 

We extracted source and background spectra for the flar- 
ing and non flaring parts for each of the four XIS units. 
The latest calibration database was used to create the XIS 
response. The response matrix files were generated using 
FTOOLS task xisrmfgen and ancillary response files were 
generated using the FTOOLS task xissimarfgen. 

To analyse the HXD data first we checked for the pro- 
cess version. Since the HXD process version was 2.0.6.13, we 
reprocessed the data as recommended by the Suzaku team. 
As a first step to reduction, we ran the command hxdtime 
to calculate HXD event arrival time correction. Then in- 
variant pulse-heights (PI) were determined using hxdpi. The 
HXD event grades were calculated using the task hxdgrade. 




2x10< 4x10' 6x10' SxlO" 

Time (in seconds) 

Figure 7. Top panel: The XISO light curve where the intensity is 
the count rate in 3-10 keV energy band. Bottom panel: overall all 
colour defined as ratio of the count rates in the 5-10 keV and 3-5 
keV energy bands as function of time. The source hardens during 
the flares. 



We used standard screening criteria (given in Suzaku ABC 
guide) to obtain the filtered PIN event files. We also ex- 
cluded the telemetry saturation period while creating fil- 
tered event files. The filtered PIN files were used to cre- 
ate the source spectra for the different flux level and dead 
time correction was performed to the PIN spectra. As sug- 
gested by the HXD team, we used version 1.2 (Method = 
LCFIT(bgd_d)) background flle to create PIN background 
spectra. Since the PIN background was simulated with a 
ten times scaled level to decrease the Poisson noise, we have 
modifled exposure keyword by entering a new exposure ten 
times the original. Since we do not have enough counts above 
30 keV, we restricted our analysis to < 30 keV and did not 
use GSO data. We added the data, background and response 
matrices from the FI (front illuminated) XIS CCDs. We used 
0.5-10 keV FI CCDs data, 0.3-10 keV BI (back illuminated) 
CCD data and 10-30 keV for PIN HXD for combined spec- 
tral fltting. We grouped the XIS data to get 3 channel per 
resolution and added 2 % systmatic error to take into ac- 
count the uncertainty in the XIS response matrix. While 
spectral fltting, we allow for calibration uncertainties be- 
tween FI CCDs, BI CCDs and PIN HXD. 

The broad band spectra of Suzaku allows for a two com- 
ponent flt to the data. We use a blackbody and Comp- 
tonization model (XSPEC model COMPTT) to represent 
the "Western" model and flt both the flaring and normal 
branch data sets. A Gaussian representing iron line emis- 
sion is statistically significant for the normal branch. Since 
the Galactic neutral hydrogen column density in the direc- 
tion of LMC X-2 is 6.3 x 10^" cm~^, we impose this as a 
lower limit to Nh- The results of spectral fits are shown in 
Table [3] The best fit Nh gets pegged to the lower limit for 
both the data sets. The spectral variation is entirely due 
changes to Comptonization. 

We repeat the spectral fitting using disk black body 
(XSPEC model DISKBB) and the same Comptonization 
model to represent "Eastern" model. The best fit spectral 
parameters are shown in Table |4l Here, the column density 
turns out to be ~ 9 x 10^° cm~^ above the galactic lower 
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Table 3. The best fit parameter obtained by fitting Suzaku data 
with black body (BBODY), Comptonization (COMPTT) and 
broad Gaussian Iron Hne models. kTgg is the blackbody tem- 
perature. kT\Y, kTg and r are the seed photon temperature, the 
electron temperature and the optical depth of the Comptoniza- 
tion model. Ep^ and ap^^ are the centroid energy and the width 
of the Iron line Gaussian. 



Para. 


NB 


FB 


Af^(102" cm"2) 


< 6.3 


< 6.3 


fcTRR (eV) 


213+? 


215+j 


kTw (eV) 


528^3 


530t* 


fcTe (eV) 


2180^6° 


2190lg" 


r 


14 05+""* 


14 75+0.04 


Ep, (keV) 


6.4 (fixed) 


6.4 (fixed) 


ap, (keV) 


'^■^'-0.14 


0.27 (fixed) 


EqWidth'' (eV) 


20 


9 


Ftest Prob': 


7 X 10-3 


0.11 


Fi 

^ line 


0.74±0.35 
463/375 


< 0.9= 
442/374 



"The column density Nu was restricted to be larger than the 
Galactic value of 6.3 X 10^° cm-^. ''The Equivalent Width of the 
Iron line. '^The F-test probability for the presence of the Gaussian 
Iron line. ''The fiux of iron line in unit of 10"^'^ erg/s/cm'^. 90% 
confidence upper limit on the iron line fiux 

limit. For the normal branch x^/d.o.f = 380/375 which is 
significantly better than the one obtained for the "West- 
ern" model X /d-o.f = 463/375. A similar trend is seen for 
the flaring branch, where the "Eastern" model gives bet- 
ter xVd-o.f = 356/373 as compared to xVd-o.f = 442/374. 
Thus spectral fitting of Suzaku data strongly prefers the 
interpretation where the soft component is disk emission 
rather than simple black body. Note that this distinction is 
possible because of the imposed lower limit on the absorp- 
tion column density. If the column density is allowed to have 
any value then a reasonable fit using the black body compo- 
nent is possible. For the model with black body emission, the 
best fit column density turns out to be Nh = 0.013 times 
10^^ cm~'^ with a x'^/dof — 367/375. For Galactic sources, 
the column density is more uncertain and hence such dis- 
tinction would not have been possible even with broad band 
data. 

The spectral fit to the non-flaring state is shown in Fig. 
[H Also plotted are the residuals with and without the addi- 
tion of a systematic error of 2%. The column density is well 
constrained to 9 x 10^° cm~^ and it is this value that has 
been used for fitting RXTE data (Table [Ijl. The inner disk 
temperature kTin, and the inner disk radius Rin as well as 
the soft photon input temperature kTw remain invariant be- 
tween the fiaring and normal branches. A broad Iron line is 
significantly detected for both branches. The electron tem- 
perature decreases during the flare, which is consistent with 
behaviour seen by RXTE, although Suzaku data requires an 



Table 4. The best fit parameter obtained by fitting Suzaku data 
with disk black body (DISKBB), Comptonization (COMPTT) 
and broad Gaussian Iron line models. kTi„ and R^bb E're the inner 
disk temperature and radius. kTyy, kTe and t are the seed photon 
temperature, the electron temperature and the optical depth of 
the Comptonization model. Ep^ and ap^, are the centroid energy 
and the width of the Iron line Gaussian. The luminosities are 
computed using the best fit model in the energy range 0.01-50 
keV and are in units of 10^* ergs/sec. 



Para. 


NB 


FB 


Nh (x 10^0 cm-2) 


8.6^?- 


8.35+»;f 


fcT,„ (eV) 


394j:^ 


3971^4 


Rdbb (km) 


99±3 


95±5 


kTw (eV) 


660 ±20 


670±20 


fcTe (eV) 


2600l^° 


2490t^° 


T 




12 40+°-l'5 


Ep, (keV) 


6.4 (fixed) 


6.4 (fixed) 


ape (keV) 


0.56l»:i| 




EqWidth" (eV) 


72 


36 


Ftest Prob*" 


7.1 X 10"^ 


2.5 X 10"* 


abs 


1.83±0.11 


1.95±0.12 




2.10±0.13 


2.23±0.12 


T ^ 

^dhb 


0.59±0.04 


0.60±0.04 


LcompTT 


1.51±0.03 


1.63±0.03 




2 60+° °^ 


2 98+''-°^ 


El"- 

Line 

xVdof 


2.68±0.54 

380/375 


1 70 + 1.42 
^- ''-'-1.04 

356/373 



"The Equivalent Width of the Iron line. ''The F-test probability 
for the presence of the Gaussian Iron line. '^Total absorbed lumi- 
nosity. ''Total unabsorbed luminosity. "^Unabsorbed luminosity of 
disk black body component. .'^Unabsorbed luminosity of Comp- 
tonized component. ^The Comptonization parameter. ''The flux 
of iron line in unit of lO"^'^ ergs/s/cm^ 



higher temperature. More significantly it is the Compton y- 
parameter which is seen to increase during the fiare, which is 
again consistent with the RXTE observations. Using a dis- 
tance of 50 kpc, the unabsorbed luminosities of flaring and 
normal branch turn out to be 2.1 x 10^* ergs/s and 2.2 x 
lO''* ergs/s, which for a I.4M0 neutron star correspond to 
1.1 and 1.15 of the Eddington value. 



4 DISCUSSION 

A detailed spectral analysis of the extragalactic source LMC 
X-2, has been undertaken using broad band data from 
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Figure 8. The top panel shows the observed Suzaku spectrum of 
LMC X-2 for the normal branch along with best fit disk blackbody 
and Comptonization model. The lower panels show the residuals 
with (bottom panel) and without (middle panel) systematic error 
addition of 2%. 



Suzaku and from a large (~ 750 ksec) data set obtained 
using RXTE-PCA. 

From the PCA data the complete Z-track evolution was 
obtained and spectra at eight representatives position were 
analysed. In the energy band of the PCA (3-20 keV) , a sin- 
gle component Comptonization fit adequately represents the 
spectra from all the track branches. This is unlik e other Z- 
sources (GX 349-H2:rAgrawal and Sreekumaij|2003l : Cyg X-2: 
iDi Salvo et"al~ll200a r where an additional soft component 
is required in the RXTE energy band. The temperature of 
Comptonizing cloud, kTe is found to decrease from the hor- 
izontal to the normal branch (s ee Fig. [H) which is sim ilar to 
what is observed for GX 17-1-2 (iDi Salvo et al. II2OO0I) but is 
in contrast to re sults from Cyg X-2 (jpi Salvo et al. |[2002l ') 
and GX 349-1-2 J Agrawal and Sreekumaijl2003l ) where the 
opposite behaviour was found. Along the Flaring branch, 
the temperature increases for LMC X-2, which is again dif- 
ferent from GX 349-1-2 where the temperature decreases. 
On the other hand, for LMC X-2 the Compton y-parameter 
decreases along the horizontal and normal branch and in- 
creases along the flaring branch. This is consistent with what 
is observed for GX 349-1-2 and it seems to be a generic fea- 
ture of Z-sources. 

These results confirm that while details of the spectral 
evolution are different for each source, there are generic be- 
haviour like the evolution of Compton y-parameter which 
seems to be similar in all Z-sources. The hard component 
evolution similarity between LMC X-2 and GX 17-1-2 sug- 
gests that there could be sub groups within Z sources where 
similar physical conditions are met. The observed similar- 
ity in the spectral evolution seen in these two sources is also 
suggested in the near identical colour and hardness intensity 
diagrams of these sources. Similar diagram are also observed 
in Sco X-1. Hence, our results are consistent with the idea 
that Z sources can be furt her divided into two subclasses , 
'Sco' and 'Cyg' like sources (| Kuulkers fc van 



I Kuulkers et al.lll997h with LMC X-2 being a member of the 
'Sco' like class. While it was suggested earlier that the dif- 
ference between ' Sco' and 'Cyg' like sou rces may be due to 
inclination angle l| Kuulkers fc van der Klia.l996i ). the spec- 
tral differences between LMC X-2 and other sources, suggest 
instead a more intrinsic difference between the two classes. 
Analysis of Z source XTE J1701-462 reveales that the source 
evolves from 'Cyg' like to 'Sco' like Z-source, suggesting 
that diffence between these two c lasses is not caused by dif- 
ference in the inclination angle l|Lin. Remillard fc HomaiTI 
l2009l : lHoman et al.ll2007l ) . 

The Suzaku data provides broad band spectral and lu- 
minosity evolution of the source from the NB to the FB. As 
expected the data requires an additional soft component. 
This additional component was represented as a black body 
(emitted from boundary layer) and as emission from mul- 
ticolour disk. The multi-colour disk approach (sometimes 
called the "Eastern" approach) provides significantly better 
fit to the data. This spectral distinction is possible because 
the Galactic neutral hydrogen column density along the di- 
rection of the source is Nh = 6.3 x 10^° cm~^ and hence the 
fitted absorption column density was restricted to be larger 
than this value. Such a constrain is usually not available for 
Galactic sources. 

The distance to LMC X-2 (50±2 kpc) is better con- 
strained than those of Galactic sources and hence broad 
band Suzaku data provides a better estimate of the bolo- 
metric luminosity. The unabsorbed luminosity of the source 
for the normal branch is 2.1 x 10^* erg/s and for the fiar- 
ing branch is ~ 2.23 x 10^* erg/s. Taking into account the 
uncertainty in the spectral fitting (~ 5%), the distance(~ 5 
%) and calibration of different instrument (< 10 %), one can 
put a highly conservative error estimate of < 15% on these 
luminosities. If the source reaches its Eddington limit on the 
normal branch, this would imply that the mass of neutron 
star is 1.50±0.1 M©. If instead the source is Eddington lim- 
ited in the fiaring branch one obtains a slightly higher mass 
of 1.61±0.09 Mq. 

The phenomenological spectral models used in this 
analysis, a single temperature Comptonized component and 
a multi-colour disk emission, is adequate considering the 
statistics of the available data and the uncertainties regard- 
ing theoretical models. However, considering that the system 
is a near Eddington neutron star source, it is expected that 
the actual radiative processes would be considerably more 
complicated with possibly quantitative and qualitative dif- 
ferences as compared to the simple model adopted here. For 
example, detailed structural analysis of the boundary layer 
reveals a temperature and d ensity stratified extended region 
(|Popham fc Sunvaev II2OO1I ). 

A fairly detailed model for Z sources is where soft pho- 
tons produced by electron cyclotron emission in the neutron 
star magneto sphere are Comptonized prim a rily in a hot cen- 
tral corona (|Psaltis Lamb fc Milled Il995l : iPsaltis fc LambI 
Il997t ). The outer accretion disk gets converted into a ra- 
dial flow which also Comptonises the outgoing photons. The 
phenomenological description in this analysis is broadly con- 
sistent with this picture with the temperature and optical 
depth of the Comptonizing component identified as the aver- 
age temperature and optical depth of the hot central corona. 
The source of seed photon for Comptonization which is as- 
sumed to be a black body with temperature, kTw in the 
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XSPEC model COMPTT, will then be the electron cy- 
clotron emission. Finally the multi-colour disk component 
is the outer disk. Note that from Table 14] the seed photon 
temperature, kTw ~ 0.65 keV is different from the inner 
disk temperature, kTin ~ 0.4 keV and hence in this inter- 
pretation the seed photons are not produced by the outer 
disk. This is consistent with the above model, where the seed 
photons are produced in the magnetosphere. 

A straightforward interpretation of the analysis is that 
the inner disk radius, Rdbb does not vary significantly when 
the source is in a flaring state and at ~ 100 km is about ten 
times the neutron star radius. The radius may be larger if 
the colour factor is significantly greater than unity as is as- 
sumed for the radius estimation. The luminosity of the disk 
emission is nearly 25% of the total luminosity and hence 
gravitational energy dissipation is not sufficient to power the 
disk. At ten times the star radius, the gravitational energy 
dissipated should be only 10% or less of the total dissipation. 
Instead the disk is probably heated by the inner corona and 
hence radiates such high luminosities. This is expected if the 
disk is truncated due to the radiation pressure of a near Ed- 
dington inner corona. However, contrary to what is seen, in 
such a scenario, one would expect the inner disk radius and 
its luminosity should vary as the source goes into a flaring 
mode. The large inferred inner disk radius argues against 
the possibility that the broad Iron line detected arises from 
the disk. Instead, it may be speculated that it arises per- 
haps from the neutron star surface. Unfortunately, the line 
is weak (with equivalent width of ~ 50 eV) and hence the 
statistics is not sufficient to test these ideas by modelling 
the line profile. 

It may also be that the phenomenological models used 
in this analysis do not sufficiently approximate the complex 
radiative process of the system. In particular, the soft seed 
photon input is assumed here to be a black body and this 
may not be the case. While the high energy spectrum is 
largely independent of the shape of the seed photon spec- 
trum, the low energy spectrum will depend on its shape. 
In the framework of a model, the magnetosphere emission 
could be more complex than a black body. Moreover there 
could be additional source of seed photons from the outer 
disk. Such complexities may vary the disk emission parame- 
ters or even eliminate the need for the component. However, 
an arbitrary seed photon spectral shape will not be falsifi- 
able by spectral analysis, since some shape will always be a 
good representation of the data. A theoretically motivated 
spectral shape for the seed photon is required. Despite these 
caveats, the spectral analysis shown in this work does reveal 
the overall behaviour of the source. Model independently one 
can still conclude that the lower energy part of the spectrum 
(which in this case is the seed photon source and the disk 
component) do not vary as the source flares and the entire 
variation can be ascribed to the Comptonizing cloud. 

The low statistics of the Suzaku data did not allow for a 
detailed broad band spectral study of the source at it evolves 
along the Z track. Instead the data was split into two broad 
parts representing the flaring and non flaring periods. This 
is unfortunate, because such a broad band evolution study 
would have been able to test whether the bolometric lumi- 
nosity increases monotonically along the track as expected 
by theoretical models and inferred by UV observations. The 
two main flares observed by Suzaku seem to be of similar 



strength (Fig[7|, however, clearly this needs to be confirmed 
by longer duration observations. This will have important 
implications on whether the source is Eddington limited dur- 
ing the flares or not. These deficiencies may be alleviated by 
a longer duration Suzaku observations undertaken simulta- 
neously with RXTE. The PCA data would distinguish the 
different track positions, while Suzaku would provide broad 
band spectra for each position. Alternatively, such results 
can be obtained by a long duration observation of the source 
by the forthcoming multi-wavelength satellite, ASTROSAT. 
The proportional counter array on board (LAXPC) would 
determine the Z-track, while the UV and other X-ray instru- 
ments would give simultaneous broad band coverage. 
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